Volume Integral Equations for Scattering from Anisotropic 

Diffraction Gratings 

Armin Lechleiter* Dinh-Liem Nguyen^ 

January 4, 2012 

Abstract 

We analyze electromagnetic scattering of TM polarized waves from a diffraction grat- 
ing consisting of a periodic, anisotropic, and possibly negative-index dielectric material. 
Such scattering problems are important for the modelization of, e.g., light propagation 
in nano-optical components and metamaterials. The periodic scattering problem can 
be reformulated as a strongly singular volume integral equation, a technique that at- 
tracts continuous interest in the engineering community, but rarely received rigorous 
theoretic treatment. In this paper we prove new (generalized) Garding inequalities in 
weighted and unweighted Sobolev spaces for the strongly singular integral equation. 
These inequalities also hold for materials for which the real part takes negative values 
inside the diffraction grating, independently of the value of the imaginary part. 



1 Introduction 

We consider scattering of time-harmonic electromagnetic waves from diffraction gratings. 
These three-dimensional dielectric structures are periodic in one spatial direction and in- 
variant in a second, orthogonal, direction (compare Figure [T]). They are used as optical 
components, e.g., to split up light into beams with different directions, and they serve in 
optical devices as, e.g., monochromators or as optical spectrometers. 




Figure 1: The diffraction grating is periodic in x\, invariant in X3 and bounded in %2- 



If the wave vector of an incident electromagnetic plane wave is chosen perpendicular to 
the invariance direction of the grating, Maxwell's equations decouple into scalar Helmholtz 
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equations, known as transverse magnetic (TM) and transverse electric (TE) modes (these 
terms are not consistently used in the literature). In this paper, we consider the equation 
of the TM mode for a non-magnetic grating, 

div (AVu) + k 2 u = 0, k > 0, 

for an a-quasi-periodic function u (that is, u(x\ + 2n,X2) = exp(27ria)u(xi, #2))). I n 
particular, we allow the real part of the discontinuous and matrix-valued material parameter 
A to be negative-definite inside the grating structure, independently of the values of the 
imaginary part. (See below for the definition of the real and imaginary - or: self- and 
non-selfadjoint - parts of a matrix.) Negative definite material parameters are a feature 
that arises in the modelization of, e.g., optical metamaterials, but also for metals at certain 
frequencies, see, e.g., [TT] . 

We reformulate the scattering problem using (a-quasi-periodic) volume integral equa- 
tions. Those turn out to be strongly singular and do not fit into the standard Riesz theory, 
since the integral operators are not compact. Nevertheless, we prove Garding inequali- 
ties for the integral equations in weighted a-quasi-periodic Sobolev spaces, which yields 
a Fredholm framework for the scattering problem. This result even holds if the real part 
Re (A) of the material parameter is negative definite inside the grating, independently of 
the imaginary part lm(A). Our approach extends a technique from |1U| . where similar 
volume integral equations have been analyzed for free space scattering problems in case 
that the scalar real-valued contrast is strictly positive. Moreover, we also prove that the 
Garding inequalities in weighted Sobolev spaces can be transformed to inequalities in stan- 
dard a-quasi-periodic Sobolev spaces, if the grating consists of isotropic material. Such 
standard Garding inequalities are particularly useful for numerical approximation since the 
corresponding Galerkin methods are easier to implement in standard Sobolev spaces than 
in weighted spaces. 

In the engineering community, volume integral equations are a popular tool to numer- 
ically solve scattering problems, see, e.g., |13}I14|. since they allow to solve problems with 
complicated material parameters via one single integral equation. However, an analysis of 
this technique in suitable function spaces is usually missing, in particular when the mate- 
rial parameters are not globally smooth, and when the arising integral operators are not 
compact. 

Recently, volume integral equations started to attract considerable interest in applied 
mathematics [5 TJ 6 | [H ]f TQ |, l 1 2 |J 1 8| . The papers [5 | ll0 p i2| analyze strongly singular integral equa- 
tions for scattering in free space. However, [12] considers media with globally continuous 
material properties, and the L 2 -theory in |10| does not yield physical solutions if the ma- 
terial parameter appearing in the highest-order coefficients are not smooth. The paper [5] 
proves a Garding inequality for a strongly singular volume integral equation arising from 
electromagnetic scattering from a (discontinuous) dielectric. Finally, [6] determines the es- 
sential spectrum of strongly singular volume integral operators arising in electromagnetic 
scattering for constant contrast. 

The analysis of the integral equation for material parameters with negative real part 
is, to the best of our knowledge, the first application of T-coercivity (a well-known frame- 
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work for variational formulations of elliptic partial differential equations with sign-changing 
coefficients, see [TH3]) to volume integral equations. As usual, the material parameter is, 
however, not allowed to take arbitrary negative values; the solvability condition for instance 
excludes that the relative material parameter takes the value —1 inside the grating. 

The paper is organized as follows: In Section [2] we briefly recall variational theory for 
the direct scattering problem. In Section [3] we derive the a-quasi-periodic volume integral 
equation. In Sections [5] and H] we prove Garding inequalities in Sobolev spaces for this 
equation. The two appendices contain two well-known results that do not fit comfortably 
into the main body of the text. 

Notation: The usual L 2 -based Sobolev and Lipschitz spaces on a domain Q are denoted 
as H s (n) and C 11 ' 1 ^), respectively. Further, #£ c (0) = {v € H S (B) for all open balls B C 
£1} . As usual, real and imaginary parts of a square matrix A are defined by Re A = 
(A + A*)/2 and ImA = (A — A*)/(2i), where A* denotes the transpose conjugate matrix. 
Both He A and Im^4 are self-adjoint and A = Re (A) + ilm(A). We denote the absolute 
value and the Euclidean norm by | • |, and the spectral matrix norm by | • |2- 

2 Problem Setting 

Propagation of time-harmonic electromagnetic waves in an inhomogeneous and isotropic 
medium without free currents is described by the time-harmonic Maxwell's equations for 
the electric and magnetic fields E and H, respectively 



where to > denotes the angular frequency e is the positive electric permittivity, no is 
the (scalar, constant and positive) magnetic permeability, and a is the conductivity. The 
permittivity and conductivity are allowed to be anisotropic, but required to be of the special 
form 



with real and symmetric 2x2 matrices et = (e«)i,i=i 2 an d o~t = j=l,2, and real func- 
tions £33 and (733. Furthermore, we assume in this paper that all three material parameters 
are independent of the third variable X3 and 27r-periodic in the first variable x\. Moreover, 
e equals £0/3 > (where I n is the n x n unit matrix) and a equals zero outside the grating. 

If an incident electromagnetic plane wave independent of the third variable X3 illu- 
minates the grating, then Maxwell's equations (pQ) for the total wave field decouple into 
two scalar partial differential equations (see, e.g., [7]). Indeed, since both, E and H 
do not depend on £3 it holds that c\xv\E = (d2E^,—diE^,diE2 — <92-Ei) T and cmlH = 
{82Hz, —d\ H3, d\H2 — c?2-f^i) T - If we plug these two relations into the Maxwell's equa- 
tions (HJ) we find that H% satisfies the two-dimensional scalar and anisotropic equation 



curl H + iujeE = aE 



curlE 1 — iuifiQH = 






m 



2 



(2) 



with wave number k := ui-^/eqJIq and material parameter 
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The usual jump conditions for the Maxwell's equations imply that the field u and the co- 
normal derivative v ■ e~ 1 S7u are continuous across interfaces with normal vector v where e r 
jumps. Note that e r is 2-7r-periodic in x\ and equals I 2 outside the grating. 

We seek for weak solutions to ([2]) and assume that e T £ L°°(IR 2 ,C 2x2 ) takes values 
in the symmetric matrices, and that e^ 1 € L°°(M?, C 2x2 ). Moreover, we suppose that 
Re(e 1 T 1 ) is pointwise strictly positive or strictly negative definite, and that Imfe^ 1 ) is 
almost everywhere positive semidefinite (even if we do never exploit the latter assumption). 
Note that we do not assume that Re (e" 1 ) is positive definite in all of M 2 . 

For the two-dimensional problem ([2|), incident electromagnetic waves reduce to u % {x) = 
exp(ifca; • d) = exp(ife(zi<ii + X20I2)) where \d\ = 1 and d 2 7^ 0. When the incident plane 
wave u l illuminates the diffraction grating there arises a scattered field u s such that the 
total field u = u l + u s satisfies ([2]). Since Au l + k 2 u % = 0, the scattered field satisfies 

div {e^Vu 3 ) + k 2 u s = -div (QVu*) in R 2 , where Q := e' 1 - I 2 (3) 

is the contrast. Note that u l is a-quasi-periodic with respect x±, that is, 

u l (x\ + 2tt,x 2 ) = e 2ma u l (xi,x 2 ) for a := kd\. 

Since u % is a-quasi-periodic and s Y is periodic, the total field and the scattered field both 
are also a-quasi-periodic in x\. For uniqueness of solution, the scattered field additionally 
has to satisfy a radiation condition. Here we require that u s above (below) the dielectric 
structure can be represented by a uniformly converging Fourier (-Ray leigh) series consisting 
of upwards (downwards) propagating or evanescent plane waves, see [4,9j, 

u s (x) = J2ufe iajX1±il * jiX2 ~ P \ *2^±P, aj:=j + a, := (k 2 - a 2 ) 1 ' 2 , (4) 

where p > sup{|x2| : (xi,X2) T G supp(Q)}. Let us set F± p = (— vr,7r) x {±p}. The numbers 
iij are the so-called Rayleigh coefficients of u s , defined by 

uf = — / u s (xi, x 2 ) exp(— ioijXi) ds , j 6 Z. 
3 271 Jr ±p 

A solution to the Helmholtz equation is called radiating if it satisfies (HJ). If k 2 > a 2 then 
the jth mode exp(iajXi ± i/3j(x 2 ± p)) is a propagating mode, whereas k 2 < a 2 means that 
exp(iajXi ± i/3j(x 2 ± p)) is an evanescent mode. 

Variational solution theory for the scattering problem ©-(JH) is well-known, see, e.g., [U 
Setting 

Q p := (-TV, n) x (-p,p) 

for p > sup{|a;2| : {x\,x 2 ) T E supp(Q)}, one can variationally reformulate the problem in 
the space H^(^l p ) := {u 6 H 1 (Q p ) : u = UW for some a-quasi-periodic U G ^^(M 2 )}. 
The resulting variational formulation is to find u s 6 H^Qp) such that 

f (e^Vu s ■ Vv-k 2 u s v)dx - f vT + (u s )ds - f vT~{u s )ds 

(5) 

QVu { ■ Vvdx 
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for all v G H^(0, p ). The operators T^, <p h-> i£) 3 - ez e^' 331 , are the so-called exte- 
rior Dirichlet-to-Neumann operators on T±„. The sesquilinear form in ([5]) is bounded on 
H^(Q p ) and satisfies a Garding inequality if, e.g., Re(e~ 1 ) is positive definite, that is, 
£*Re (e~ )(#)£ > c|£| 2 > for £ G C 2 and almost every x G In this case, analytic 
Fredholm theory implies that the set of real wave numbers (excluding Rayleigh frequen- 
cies) where non-uniqueness occurs is at most countable, see [HE]. If Re (e^ 1 ) changes sign, 
proving Fredholm properties of the variational formulation ([5]) is non-trivial, at least if 
Ime" 1 vanishes. 

In this paper, we establish a Fredholm framework for the scattering problem via integral 
equation techniques, that is, uniqueness of solution implies existence. We do not aim to 
prove the corresponding uniqueness results, since for periodic scattering problems such 
results are anyway not available at all frequencies, except under restrictive geometric (non- 
trapping) conditions, see [I]. 



3 Integral Equation Formulation 

In this section, we reformulate the scattering problem ([2]) as a volume integral equation, 
and prove mapping properties of the integral operator between Sobolev spaces. To this end, 
let us recall that Q = e~ l — I2 and denote by D C Q p the support of Q, restricted to one 
period {— tt < x\ < it}. By Gk, a we denote the Green's function to the a-quasi-periodic 
Helmholtz equation in M 2 , see [9j. Under the assumption that 

k 2 ^ a 2 for all j G Z, (6) 

the a-quasi-periodic Green's function has the series representation 

G k , a (x) ■= ^2 ^-exp(ia j xi + i/3 i |x 2 |), x = G M 2 , x / {^^j for meZ. (7) 

Due to © all the /3j = (k 2 — a 2 ) 1 / 2 are non-zero. 

Lemma 1 (Cf. [9]). The Green's function G^^ a can be split into Gfc jQ (x) = {1/ 4)H^ (k\x\) + 
ty(x) in M. 2 where \l/ is an analytic function solving the Helmholtz equation + k 2 ^ = 
in (-2vr,2vr) x R. 

We also define a periodized Green's function, firstly setting 

IC p (x) := G k , a (x), x = G R x (-p,p), x ^ ( 2 ™J for m G Z, (8) 

and secondly extending K p {x) 2p-periodically in X2 to M 2 . 
The trigonometric polynomials 

1 (•(■ . x . -hit ■ \ ■ fji\ ^2 



¥j(x) ■= exp + a)xi + 1—^2 J , J ' = ( ^ ) G Z , (9) 
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are orthonormal in L 2 (£l p ). They differ from the usual Fourier basis (see, e.g., |16[ Section 
10.5.2]) only by a phase factor exp(iaxi), and hence also form a basis of L 2 (Q p ). For 
/ G L 2 (Sl p ) and j = (ji, J2) T £ %> 2 , f(j) '■= Jq f^Pjdx are the Fourier coefficients of /. 
For < s < oo we define a fractional Sobolev space Hp er (Q p ) as the subspace of functions 
in L 2 (Q p ) such that 

ii/ii^ p) = E( 1 + ui 2 ) s i^)i 2 < 00 - ( 10 ) 

It is well-known that for integer values of s, these spaces correspond to spaces of a-quasi- 
periodic functions that are s times weakly differentiable, and that the above norm is then 
equivalent to the usual integral norms. 

Lemma [T] implies in particular that tC p has an integrable singularity and that the Fourier 
coefficients IC p (j) are well-defined. To compute these coefficients explicitly, we set 



A,: lr (j\ ha) 2 tor,G^ 



Theorem 2. Assume that k 2 ^ a 2 for all j G Z. Then the Fourier coefficients of the kernel 
K, p from (JHJ) are given by 



Pup 



cosQ^TrJe 1 J l 
■\ZAttp Xj 
i_ /P\ 3 / 2 

4j2 V 7T / 



for Xj 0, 
else, 



Remark 3. Note that K, p (j) is well-defined for Xj = 0: Since k 2 ^ a 2 for all n G 7L, the 
definition of Xj implies that ji 7^ whenever Xj = 0. For completeness, we include a proof, 
noting that the case Xj ^ is also shown in [lSi Section 1.1]. 

Proof. It is easy to check that (A + k 2 )tpj = Xjifj for j = (ji,j2) G 1? . If Xj ^ 0, Green's 
second identity implies that 

Mi) = 



JC p (x)ipj(x) dx =X 



AT 1 lim 




Gk,a{ x ){^- + k 2 )(fj(x) dx 



■<P3 



ds 



Q P \B(0,8) 



(A + k )Gk, a (x) ( Pj{%) dx 



(11) 

(12) 



where v denotes the exterior normal vector to B(0,5). The last volume integral vanishes 
since (A + k 2 )Gk, a = in Q p \ -B (0, 6) for any 5 > 0. Let us now consider the first integral 
in (jlip . The boundary of £l p consists of two horizontal lines T± p and two vertical lines 
{(xi,X2) : xi = ±7r, — p < X2 < p}. Hence, the normal vector v on these boundaries is 



either (±1,0) T or (0,±1) T . Straightforward computations yield that 

G k , a (x lt ±p) = — e -"*i, 9 2 G k , a ( Xl ,±p) = T — ^e^" 11 , (13) 

1 _■ ij2'T 

^•(xi,±/>) = -^=e ^i" 1 cos(j 2 7r), and d2<Pj(xi,±p) = —(p(xi,±p). (14) 

In consequence, 

/ \Gk,or^- - dG ^' a j Pj \ ds = - [ d 2 G k , a Tp]ds + f d 2 G k , a 7p]ds 
Jr ±p \ ov Ov J j Tp j r _ p 

= -2 / d 2 G kta if]ds . 
Using the above formulas for d 2 G kj0t and TpJ in (|13|) and (JHJ), respectively, we find that 

Computing the partial derivatives of G k , a and ipj with respect to x\ analogously to the 
above computations, one finds that the integrals on the vertical boundaries of Q p vanish 
due to the a-quasi-periodicity of both functions. Thus, we obtain that 

f ( n dG k,a—\ , cos(j 2 vr) i/3 

Now we consider the second integral in (jlip . From Lemma [T] we know that G k a (x) = 
jHq (k\x\) + ty(x) where ^ is a smooth function in fl p . Obviously, 

, im / (*M » w ) ds=0 . 

The asymptotics of and its derivative for small arguments, 

HPM = -logr + 0(l) and (HP)' (r) = — + 0{1) as r 0, 
allow to show that 

]i- / - ^N?) d-,= -^=, (16) 



see, e.g., |16t Theorem 2.2.1]. Combining (|15p with (|16|) yields that 

K p {j) = -^-(cosih^nP - l) for \j ± 0. 
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For Xj = we use de L'Hospital's rule to find that 

cos(j 2 7r) exp(ipy/j- (ji + a) 2 ) - 1 _ ip 3 / 2 



7^(ii+«) 2 +0 2 7r/p) 2 V^Trp [7 - (Ji + a) 2 - IWp) 2 ] 4vr 3 / 2 j 2 ' 

Note that the assumption that k 2 ^ a 2 for all j E Z 2 implies that Xj and j 2 cannot vanish 
simultaneously. □ 

Since the Fourier coefficients of the /C p decay quadratically, \K- P (j)\ < C/(l + (ji + a) 2 + 
(j27r/p) 2 ) for j E Z 2 , the convolution operator with kernel JC p is bounded from L 2 (Q p ) into 

-^per(^p)- 

Proposition 4. Assume that k 2 ^ a 2 for all j E Z. T/ien i/ie convolution operator K p , 
defined by 

(K p f)(x)= K p {x-y)f{y)dy for x E Q p , 
Jn p 

is bounded from L 2 (Q, p ) into H 2 CT (Q p ). 

Recall that D C £l p is the support of contrast Q. Let us additionally introduce 

n ■= (-Ti, -k) x r 

and H^(VIr) := {u E H e (tt{i) : u = for some a-quasi-periodic C7 E //^(R 2 )} for 

^ E N, i? > 0. 

Lemma 5. Assume that k 2 7^ a 2 /or aZZ j E Z. T/ien f/ie volume potential V k defined by 

(Ykf)(x)= / G kfCt (x-y)f(y)dy, x £ fl, 
Jd 

is bounded from L 2 (D) into H 2 (Qr) for all R > 0. 

Proof. Consider x £ C°°(T2) such that x = 1 i n D, < x < 1 in fi p \ D and x( x ) = for 
\x 2 \ > p. Then V k g = X V k g + (1 - x)^. Note that (1 - x )V k g = f D (l - X )G(- - y)g{y) dy 
is an integral operator with a smooth kernel, since the series in ([7]) converges absolutely 
and uniformly for \x2\ > p > 0, as well as all its partial derivatives. In consequence, the 
integral operator (I — \)Vk is bounded from L 2 (D) into H 2 (Qji), since 



\d^d^((l- X )V k g)\\ 2 L2{QR) < [ [ \d^d^[{l- X (x))G k , a (x-y)]\ 2 dydx \\gf„ 

y ' J flu J D 



(D) 



for all /3i t 2 E N such that f}\ + /3 2 < 2. It remains to show the boundedness of x^4 from 
L 2 (D) into H 2 (Q, p ). Let (7 E L 2 (D) and consider the operator ET 2p from Proposition 21 
mapping L 2 (Sl 2p ) into H 2 CT (n 2p ) C F 2 (0 2p ), 

(K 2p g){x) = j JC 2p {x - y)g{y) dy for x E fi 2 p- 

If x E ftp, then |x 2 — y 2 | < 2p, that is, /C 2p (x — y) = G k ^ a (x — y). Hence, K2 P g = V k g in 
Q p , and hence X K2 P g = xVkg in ^p- Since x is a smooth function, we conclude that xVk is 
bounded from L 2 (D) into #„(fi p ). □ 
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Note that the potential V k f can be extended to an a-quasi-periodic function in H? 
due to the a-quasi-periodicity of the kernel. 

Lemma 6. For g £ L 2 (D,C 2 ) the potential w = divV^g belongs to H^(Q p ) for all p > 0. 
It is the unique radiating weak solution to Aw + k 2 w = —divg in Q,, that is, it satisfies 

/ (Vu) • VU — k 2 wv) dx = — I g-Vvdx (17) 
Jn Jd 

for all v £ H^(Q) with compact support, and additionally the Rayleigh expansion condi- 
tion dU). 

Proof. Lemma [5] and a-quasi-periodicity of the kernel of V k imply that w is a function in 
Ha(Q p ) for all p > 0. It is sufficient to prove (fT7]) for all smooth a-quasi-periodic test 
functions v that are supported in {\x2\ < C} for some C > depending on v. It is 
well-known that p = Vug £ H 2 (Q) is a weak solution to the Helmholtz equation, that is, 



(Vpj • VdjV — k 2 pjdjv) dx = — gjdjvdx 
i ' ' Jd 

for j = 1,2. An integration by parts shows that 

/ (Vdivp • X7v — k 2 d'wpv) dx = — I g-Vvdx, 
Jn Jd 

which implies (fT7|) due to divp = w. Since the components of the potential p = V^g satisfy 
the Rayleigh condition, a simple computation shows that the divergence w = divp does 
also satisfy the latter condition. Uniqueness of a radiating solution to (|17p when g = 
follows from the Rayleigh expansion condition and a unique continuation argument. □ 

Returning to the differential equation Q for the scattered field u s , let us set / = QX7u l . 
(Recall that Q = e" 1 — J2O The variational formulation of © is 

I (Vm s • VU - k 2 u s v) dx = - [ (QVu s + /) • Vvdx (18) 
Jn Jd 

for all v £ H^i^l) with compact support in Q. From Lemma [6] we know that the radiating 
solution to this problem is given by u s = divVfc(QVn s + /). Hence, we aim to find u s : 
O — > C that belongs to H^(Qr) for all R > 0, such that 

u s - dwV k (QVu s ) = dwV k (f) in fl (19) 



4 Garding Inequalities in Weighted Sobolev Spaces 

For scattering problems in free space and for scalar and positive contrast, the paper |10j 
investigates integral equations similar to (|19p in weighted spaces. In this section we general- 
ize the results from [10] to anisotropic and possibly sign-changing coefficients in a periodic 
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setting, proving a Garding inequality for / — divVfc(QV-) in a anisotropically weighted 
a-quasi-periodic iJ 1 -space. 

From (|19p it is obvious that the knowledge of u in D is sufficient to determine u in 
Q \ D by integration. Thus, we define the operator : f \— > div Vkf that is bounded from 
L 2 (Z?,C 2 ) into H^(D) and consider the integral equation 

u = L k {QVu + f) in (20) 

To study Garding inequalities for volume integral equations, we introduce suitable weighted 
Sobolev spaces. To this end, we recall that the symmetric 2x2 matrix Re (Q) has pointwise 
almost everywhere in D an eigenvalue decomposition Re (Q) = U*YiU with a diagonal 
matrix £ and an orthogonal matrix U. This decomposition can be used to define the 
absolute value |Re(Q)| = U*\T\U and the square root |Re {Q)\ 1/2 = C/|X| 1/2 £7*, where the 
absolute value and the square root are element-wise applied to the diagonal matrix S. The 
two eigenvalues Ai j2 of Re (Q) define 

Amin(^) = min{|Ai(x)|, |A 2 (x)|}, A max (x) = max{|Ai(x)|, |A 2 (x)|}, x £ D. (21) 

We assume in the following that Re (Q) is pointwise either strictly positive or strictly 
negative definite, such that we can assign a sign function sign(Re (Q)) G L°°(0) to Re (Q), 
indicating whether the eigenvalues of Re (Q) are positive or negative. In the sequel, we 
write Re (Q) > c in D (Re (Q) < c in D) to indicate that the eigenvalues Ai j2 are larger 
than or equal to (less than or equal to) a constant c, almost everywhere in D. Note that 
the spectral matrix norm is denoted by | • | 2 . 

We denote by q{D) the completion of H^(D) with respect to the norm || • \\ H i q (d)-> 

\HhI iQ (d) ■= W\Xe(Q)\Vv\\h(D,&) + Hh(Dy ( 22 ) 

Since we assumed that supp(Re(Q)) = D, this norm is non-degenerate. Moreover, 
is an equivalent norm in H^(D) provided that |Re (Q)\ is bounded from below in 

D by some positive constant. In general, ^ < (l+|||-\/|Re (Q)|| 2 ||z°°(d)) IMI-ffi(-D)- 

Note also that the norm of q{D) is linked to the sesquilinear form 

a Q ( u ,v) = [sign(Re (Q))QVu ■ Vv + uv] dx , u,veH^ Q (D). (23) 
Jd 

Indeed, ||u||^i (D) = Re [clq(u,u)] for u £ H^q(D). In consequence, the form aq is 

non-degenerate, that is, if ag(u, v ) = for all v 6 n(D), then u = 0. 

If Im.Q vanishes in D (that is, the values of x i— > Q(x) are self-adjoint matrices), then 
clq is simply the inner product associated with the norm of H\ q(D), 

{u,v) h i q{d) = [|Q|V«- Vv + uv] dx, u,veH\ >Q {D). 
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Lemma 7. Assume that there exists C > such that 

|Im (Q(ac))C| < C|Re (Q(x))f | /or a/mosi euen/ x e D and all (eC 2 . (24) 
T/ien u i— > Lk(QS7v) is bounded on q(D). 

Proof. Due to Theorem [5l L k is bounded from L 2 (D,C 2 ) into H^(D). Furthermore, v i— )■ 
QVw is bounded from H^q(D) into L 2 (D,C 2 ), since 

IIQVxt|| jL 2 (jDjC 2 ) < ||Re (Q)Vii|| i 2( £ j > c 2 ) + ||Im (Q)Vii|| jL 2 (ZP:C 2 ) 

< |||Re(Q)|V«|| £ ,2 (AC 2 ) +C||Re(Q)V«|| i 2 (DiC 2) (25) 

< (1 + C) || I VlRe (Q)| \ 2 \\l<~P) \Hh^ q (d) ■ 

Moreover, the imbedding H^(D) C H\q{D) is bounded, as mentioned above. Hence, 
v 1 y LkiQVv) is bounded on n(D). □ 

Remark 8. Condition (|24p is satisfied if the absolute values of the eigenvalues oflmQ 
are pointwise bounded by CA m i n (recall from ([2T]) f/mf A m i n zs i/ie minimum of the absolute 
values of the eigenvalues of Re (Q) ). 

If u G H\{D) C H^q(D) solves the Lippmann-Schwinger equation (f20j) . then 
Lemma (J7J) implies that u solves the same equation in n(D). Since clq is non-degenerate, 
solving the Lippmann-Schwinger equation in H^q{D) is equivalent to solve 

a Q (u-L k (QVu + f),v) = for all v G H^ q(D). (26) 

If u G H^q(D) solves the latter variational problem for some / G L 2 (D,C 2 ), then u = 
L k (QVu + /) belongs to #£(£>), due to (J2SJ) and since is bounded from L 2 (D, C 2 ) into 

Proposition 9. Assume that f G L 2 (D,C 2 ). Then any solution to the Lippmann- 
Schwinger equation (|2U|) in H^(D) is a solution in H^q(D) and vice versa. 

Our aim is now to prove a (generalized) Garding inequality for the variational prob- 
lem (|26p . The following lemma will turn out to be useful. 

Lemma 10. Suppose that X and Y are Hilbert spaces. Let T\ t 2 be bounded linear operators 
from X into Y and consider the sesquilinear form a : X x X — > C, defined by a(u, v) = 
{Tiu,T2v)y for u,v G X. If one of the operators T\ and T2 is compact, then the linear 
operator A: X — > X, defined by {Au,v)x = a(u,v) for all u, v G X, is compact, too. 

Proof. It is easily seen that A is a well-defined bounded linear operator. Obviously, 
\(Au, v)x\ = \a(u, v)\ < C||Tiit||y||T2v||y for u,v G X. Assume that T\ is compact, and 
note that 

II A II \{Au,V)x\ . ^urr I, 

\\Au\\x = sup — 7— 7 < C||Tin||y. 

o^v&x \\v\\x 
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If a sequence {u n } converges weakly to zero in X, then {Tiu n } contains a strongly con- 
vergent subsequence tending to zero in Y. Consequently, {Au n } also contains a strongly 
convergent zero sequence, which means that A is compact. One can analogously derive the 
compactness of T in case that T2 is compact, since a(u,v) = (T^TiUjV). □ 

The next lemma proves Garding inequalities for the operator - L^QVv) using 

the sesquilinear form oq from (|23p . The second part of the claim uses a periodic extension 
operator 

E : Hl(D) -> Hl(n), E{u)\ D = u, E{u)\ n ^ 2p = 0, 
introduced in Appendix [Bj The operator norm of E is 

\\E\\m a (D)^m a (n 2p ) = (l + \\ E \\ 2 H i { D)^m a (n 2p \D)) ■ 

Theorem 11. Assume that D is a Lipschitz domain and that Q G L°°(D, C 2x2 ). 

(a,) IfKe(Q) > in D, then there exists a compact operator K + on H^q(D) such that 

Re [a Q (v - L k (QVv), v)} > \\v\\ 2 h ^ q{d) - Re {K+v, v) h x q{d) , v G ^(D). (27) 

(b) If Re (Q) < -1, and j/ 

l^[|Hi(i?HflJ(n 9p ) < mf |Re(Q)|f , (28) 

i/ien t/iere exists a constant C > and a compact operator K- on H^q(D) such that 

-Re [a Q (v - L k (QVv), v)} > C\\vf Hl , D) -Be (K.v, v) Hl {D) , v G H^ Q (D). (29) 

Remark 12. // Im (Q) = in D, then both statements (|27p and (|29p are nothing but 
standard Garding estimates: The form oq defines an inner product on H^q(D), and, 
e.g., (|27|) can 6e rewritten as Re — L k (QVv), > ||w|| 2 — Re {K + v, v) for v G id^ q(D). 

Proof, (a) We start with the case Re (Q) > in D. Let v G id* q(-D) and define w by 

to = UiQVv) = div / Gi, a (--y)[Q(y)Vv(y)]dy in fi. (30) 

Then w G id*(Q) decays exponentially to zero as |x2| tends to infinity. Moreover, Aw — w = 
—div (QVu) holds in O in the weak sense due to Lemma El that is, 

/ [Vio • + imp\ dx = - / QVvV^dx for all ^ G id^(fi). (31) 

Setting V> = to, we find that —Re f D QS/v ■ S/wdx = IMI^i^- Hence, 

Re [a Q (v- Li(QVv),v)]= / [Re (Q)Vv ■ VU + \v\ 2 ] dx - Re / [QVw • VU + wv] dx 

Jd Jd 

= / [\y / Re{Q)Vv\ 2 + \v\ 2 -Re{wv)]dx + / [|Vn;| 2 + \w\ 2 ] dx 
Jd Jn 

1 1 f 

>M\ 2 H i aQ{D) --\\v\\ 2 L 2 {D) + - J^[\v\ 2 + \w\ 2 -2Re(wv)]dx , 
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where the last term on the right is positive. In consequence, 

Re [a Q (v - L k (QVv), v)] > \\v\\ 2 Hl (D) - ~(v, v) L 2 (D) - Re [a Q ((L k - L{)(QVv), v)] 

for all v £ n{D). Due to Lemma [TU1 and Rellich's lemma there exists a compact operator 
K\ on H\q{D) such that (v,v) L 2^ = 2Ke(Kiv, v) H i q (d)- Further, the operator (L k — 

Li)(QV-) is compact on H^(D) due to the smoothness of the kernel shown in Appendix [XI 
Hence the operator K2 defined by (K2V, v) H i (D \ = an((Lk — Li)(QVv), v) is compact 

on H^q(D) due to Lemma flQl and the boundedness of the imbedding H^(D) C H\q{D). 
Setting K + := K\ + if 2, we obtain the claimed generalized Garding inequality. 

(b) Now we consider the case that Re (Q) < —1 in I?, and assume additionally that (|28p 
holds. As in the first part of the proof, the variational formulation (|3ip for w, defined as 
in (f30|) . yields that 



-Re [a Q (v - Li(QVv), v)} = Re / [Re (Q)Vv ■ VU - \v\ 2 - QVw ■ Vv + wv] dx 

Jd 

= -j [| v / |Re(Q)|Vu| 2 + \v\ 2 ] dx + |MI#i(n) +Re J wvdx 



> IIHIfil(n) ~ II U H^ Q (D) + Re j D ^vdx. 



We plug in ip = —E{v) into pip and take the real part of that equation, to find that 
|Re(Q)|Vw||| 2(DjC2) < |MlHi(n)ll-E(w)llffi(n) < \\E\\ H i {D) ^ H i {n2p) \\w\\ H i {n) \\v\\ H i {D) 



--1 



< ll-^ll lk||fli(n) y||V|Re(Q)| \2\\l°°(d)\MhI q (d) + IMIl2(d) 

For x £ D, the spectral matrix norm | y^Re (Q)\ (x)\2 of the inverse of \/|Re (Q)|(x) 
equals the reciprocal value Ami^x) -1 / 2 (A m in imax are the smallest /largest eigenvalue, in 
magnitude, of Re(Q), see (f2Tj) ). Note that 

\\\V\^W\~ l \2\\li( D ) = [supA^n^)- 1 / 2 ]- 1 = inf A^nOr) 1 / 2 < SU P \^(x) 1/2 
( ' x£D X <^ D xGD 

< [l + su P A max (x)] 1 / 2 = [l + |||Re(Q)| 2 || Loo(D) ] 1 / 2 . 

Next, we estimate that 

^h, Q (D)~ [ 1+ \\\ Re (QML^(D)]\\v\\ 2 L 2 {D) < \\v\\ 2 h 1 q(d) -\\v\\ 2 L 2 {d) 



< \\ E \\ WW |Re(Q)l |2|U°°(z>) \\w\\h^) 

(\\ v \\hI q (D) + [1 + ll|R-e(Q)|2||L°°(D)] 1/2 |bllL2(D)) • 

Dividing by the term in brackets on the right, we obtain that 
WvWh 1 a (D)~ [ 1 + \\\^(Q)\2\\oo] 1/2 \\v\\ L 2 {D) < ||£|| |||v^|Re(Q)| \ 2 \\l°°(d) ||Hlfl*(n)- ( 32 ) 
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Note that the constant 



c : ~ \\ E \\Hl(D)-+Hi(n 2p ) \\W\R e (Q)\ h\\L°°(D) 

is by assumption (f28|) less than one. If we set for a moment, C = [1 + |||Re (Q^Hoo] 1 / 2 
then ([32]) and Cauchy's inequality imply that 

c2 W w \\ 2 Hi(n) > \M 2 H i aQ{D) + C 2 \\v\\ 2 L2{D) -2C\\v\\ H i aQiD) \\v\\ L 2 iD) 

> (1 - e 2 )\\v\\ 2 Hi Q(D) + C 2 (l - l/e 2 )\\v\\ 2 L2{D) , e G (0, 1). 

In consequence, 

- Re [a Q (v - L k {QX7v), v)] > ~ lj \\vf H i q{d) 

f e 2 - 1 

-Re / wvdx +C 2 —- T \\v\\ 2 L2(D) + Re[a Q ((L k -L i ){QVv),v)} (33) 

Jd \ ce ) 

for e G (0, 1). Since c < 1 there exists e G (0, 1) such that 1— e 2 > c 2 , that is, (1— e 2 )/c 2 — 1 > 
0. The last three terms on the right-hand side of (|33p can then be treated as compact 
perturbations, in a similar way as in the proof of the first part. □ 

Remark 13. (a) IfKe(Q) < —1 in D, then solutions to div((/2 + Q)Vu) + k 2 u decay 
exponentially in D. If not only the electric permittivity but also the magnetic permeabil- 
ity changed sign, then the corresponding solution would not decay, yielding a possibly more 
interesting metamaterial. Volume integral equations for such structures yield operator equa- 
tions combining L k and V k , see, e.g., JiflJ /. Since V k is compact on H , the above Gdrding 
inequalities extend to this setting. For simplicity, we restrict ourselves here to the non- 
magnetic case. 

(b) In the last result, we assumed that the sign of He (Q) is constant in D. It is possible 
to treat sign changes of the contrast function in D, but the simple choice ip = —E(v) that 
we plugged in the second part of the proof into (|3ip has to be adapted. 

It is a standard result that the Garding inequalities from the last theorem imply the 
following consequences for the solvability of the integral equation and the scattering prob- 
lem. 

Theorem 14. Suppose that the assumptions of Theorem \llVa) or (b) hold, that the 
boundedness condition (f24"|) holds, and that the homogeneous equation v — L k (QVv) = 
in H^q(D) has only the trivial solution. Then (|20p has a unique solution for all 
f G L 2 (D,C 2 ). If f = QVu\ then this solution can be extended by the right-hand side 
of (|20p to a solution to the variational formulation of the scattering problem ([5]). Espe- 
cially, if the integral equation is uniquely solvable in H^q(D), then ([5]) is uniquely solvable 
inH l a {VL p ). 
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5 Garding Inequalities in Standard Sobolev Spaces 



The generalized Garding inequalities from the last section imply Gfarding inequalities in 
the standard unweighted periodic Sobolev space H^(D) if the material parameter e r (or, 
equivalently, the contrast), is isotropic. Hence, in this section we assume that the contrast 
is a scalar real-valued function q, that is, 

Q = ql2 in fi. 

As above, D is the support of q. Under this assumption we denote the weighted Sobolev 
spaces from ([22]) by H^ q (D), and their norm by 

IHh^D) ■= (lM Re (tf)|Vw||£ 2(AC2) + IMlia (J ,)) ' ■ 

Since q is real- valued, the form a q from ([23]) is the inner product of H^JD), and the 
generalized Garding inequalities from the last section directly transform to standard ones. 
Again, we assume that the sign of q is constant in D. Since we use regularity theory to 
prove compactness of certain commutators, we will need to require more smoothness of q 
and D compared to the results in the last section. 

Lemma 15. Assume that D is a domain of class C 2,1 and that [i G C 2,1 (D), 2ir -periodic 
in x\. Then T : H^(D) — > H^(D) defined by Tv := div [fiVk(qV (v / /j,)) — Vk(qVv)] is a 
compact operator. 

Proof. We denote by /x* G C 2,1 (f2 p ) a periodic extension of \i G C 2,1 {D) to Q p (see Ap- 
pendix [B] on periodic extension operators). Then /j,*\d = M- Consider the two a-quasi- 
periodic functions 

w\ = Vk(qV(v/fi)) and w 2 = Vk{qVv) in tt p . 
Both functions satisfy differential equations, 

A^ Wl ) + k 2 ^ Wl ) = l- q ^ {v/ ^ + 2V ^ Vwi + WlA ^ inA _ 

[2Vn*- Vwi + wiAfi* in Sl p \ L>, 

and A1U2 + k 2 u>2 = —qVv in D and Au^ + k 2 W2 = in il p \ D. Hence, w = ^*w\ — W2 
solves 

a ,72 J-g/iV(l//i)u + 2V/x- Vwi + wxA/j, =: gx in D, 

Aw + k w = < 

\wiAn* + 2X7 fi* • Vu>i =: g 2 in Q, p \ D. 

The functions g± and g2 belong to H^(D) and H^(Q p \ D), respectively. Their norms in 
these spaces are bounded by the norm of [x in C 2,1 (D) times the norm of v in H^(D). Due 
to LemmaO the jump of the trace and the normal trace of wi } 2 across dD vanishes. Hence, 
the Cauchy data of w are also continuous across the boundary of D. 

Since the volume potential is bounded from L 2 (D) into H 2 (D), it is clear that w 
belongs to H 2 (D). The smoothness assumptions on D and fi moreover allow to apply elliptic 
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transmission regularity results Theorem 4.20] to conclude that w is even smoother 
than H 2 . These regularity results will in turn imply the compactness of the operator 
T : v i— > divw on H^(D). A straightforward adaption of the transmission regularity 
result Theorem 4.20] to the periodic setting shows that 



\ W \\H 3 (D) 



< C 



ilJ 



\w(n p ) + \\gi\\m(D) + \\92\\ H i {np \D) 



< C\\v\ 



HUD)' 



□ 



The following lemma shows that the Garding inequalities in the weighted spaces H^ q (D) 
can be transformed into estimates in H^(D) if, roughly speaking, the real- valued contrast 
q is smooth enough and if (Vq)/q is bounded. 

Theorem 16. Assume that the scalar contrast q is real-valued, that \q\ > q$ > in D, and 
that y/\q\ £ C 2 ' l (D). Moreover, assume that D is of class C 2 ' 1 . 

(a) If q > there exists a compact operator K + on H^(D) such that 



Re (v - L k (qVv), v) H i {D) > \\v\\ 2 H i {D) - Re {K+v, v) H i {D) , 
(b) lfq<0, and if 



v e Hl{D). 



(34) 



\\E\\m a (D)^m a (n 2p ) < inf \q\ 1/2 , 
then there exists a compact operator K- on H^(D) such that 

-Re{v - L k (qVv),v) H i {D) > C\\v\\ 2 H i {D) - Re{K_v, v) H i {D) , v G H*(D), 
where C is the constant from ([291) . 

Proof. We only prove case (a) here, supposing that q > qo > in D. The proof for case (b) 
is analogous, essentially one needs to replace y/q by \/\q\. For simplicity, let us from now 
on abbreviate 

fi-^qeC^iD). 

Choose an arbitrary u £ H^(D) and consider v = u/fJ,. Our assumptions on q imply 
that v E H^ q (D), since \\vf H i jD) < (2 + + 2||(Vm)/A)N&i (d) . In Th eo- 

rem fTlT a) (see also Remark [T2"]) we showed that Re (v — L k (qX7v), v) H i ^^ D j > \\v\\ 2 H1 ^ — 
Ke(Kiv, v) H i i D \ for a compact operator K\ on AD). This implies that 



Re(«- L k (qVu), u) H i (D) > \\u\\ 2 H i( D) + Re (K^u/n), u/ii) H i q{D) 



+ Re(K 2 u, V«) L 2 (AC 2) + Re (K 3 u, u) i2(D) , 



with operators 



K 2 u = V 



div [fj,V k (qV(u/fjL))-V k (qVu)] -V V//T fc (gV(u//i)) +(Vfj,)L k (qV(u/fj,)) 
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and 

K 3 u = gV(l//x) • [VL k (qV(u/fi))] + L k (qV{u/n))/(i - L k (qVu). 

Lemma [T5l the smoothness of q, and the boundedness of V k and L k from L 2 {D) and 
L 2 (D,C 2 ) into H 2 (D) and H^(D), respectively, show that and K3 are compact and 
bounded from H^{D) into L 2 (D), respectively. Then the compact embedding H^(D) C 
L 2 (D) and an application of Lemma [10] imply the claim. □ 

Remark 17. The regularity assumptions on dD and q can be weakened using more so- 
phisticated regularity results. It is for instance possible to treat piecewise smooth q by an 
analogous technique. We do not discuss this issue to avoid technicalities that would not add 
new ideas. 

Theorem 18. Suppose that q and D satisfy the assumptions of Theorem MbY a) or (b), and 
that the homogeneous equation v — L k (qVv) = in H^JD) has only the trivial solution. 
Then (|20p has a unique solution for all f G L 2 (D,C 2 ). If f = qVu 1 , then this solution 
can be extended by the right-hand side of (|20p to a solution to the scattering problem ([5]). 
Especially, if the integral equation is uniquely solvable, then (|20p is also uniquely solvable. 

A Smoothness of the Difference of Periodic Green's Func- 
tions 

The following lemma is a consequence of the corresponding result for the fundamental 
solution to the Helmholtz equation in free-space. 

Lemma 19. Assume that k 2 7^ a 2 for all j G Z. Then the difference G kja — G\^ a can be 
written as 

G M (x) - G i)Q (x) = a(\x\ 2 ) + C|x| 2 ln(|x|)/3(|x| 2 ) 
where a and j3 are analytic functions and C is a constant. 

The smoothness of the difference Gk )a —G^ a implies the following compactness statement 
for the corresponding volume potentials. 

Corollary 20. Assume that k 2 7^ a 2 for all j G Z. Then L k — L\ is compact on H^(D). 

B Periodic Extension Operators 

In this section, we exemplary show how to construct a periodic extension operator 

E : HiiD) -+ Hl(n), E{u)\ D = u, E(u)\^ n2p = 0, 

that is used in Theorem [TTJ We will only construct E for the case that the boundary 
of D = {{xi,X2) Y : x\ G (— vr,7r), C-O^l) < x 2 < (+{ x i)} is given by two 2-7r-periodic 
Lipschitz continuous functions £± : M — > (— p, p) such that £_ < — 2p/3, C+ > 2/9/3, and 
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IC±(*i) - C±K)I < M \ x i ~ x 'i\ for ^i,xi g 

patches as in [11, Appendix A]. 
For u G H\{D), we define 



v(xi,x 2 ) 



'-u(xi,2C+(xi) 
ti(xi,x 2 ) 
,w(xi,2C-(xi) 



J 



^2 J 



The general case can be tackled using local 



if C+(xi) < x 2 < 2C+(xi) - C-(xi), 

if C-(«l) < x 2 < C+(zi)> 

if 2C-(xi) - C + (iCi) < x 2 < C-(^i)- 



Note that 2£ + (xi) — C-( x i) > 2p and that 2£_(xi) — C+( x i) < —2/?. Straightforward 
computations show that |M|#i(n 2p ) < max(\/3, 2\/2M)||ti||#i fjj\, and the definition of v 
implies that this function is a-quasi-periodic. 

To define the periodic extension operator, we use a smooth cut-off function x '■ K ~~ * ^ 
that satisfies < x < 1> x( x 2) = 1 for |x 2 | < p, and = for |x 2 | > 2p. Then we set 



E(u) = w, w(x) 




for x G J7 2p , 
else. 
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